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Abstract 



Recently we have studied the Bloch effective Hamiltonian approach to 
bound states in 2+1 dimensional gauge theories. Numerical calculations were 
carried out to investigate the vanishing energy denominator problem. In this 
work we study similarity renormalization approach to the same problem. By 
performing analytical calculations with a step function form for the similarity 
factor, we show that in addition to curing the vanishing energy denominator 
problem, similarity approach generates linear confining interaction for large 
transverse separations. However, for large longitudinal separations, the gen- 
erated interaction grows only as the square root of the longitudinal separation 
and hence produces violations of rotational symmetry in the spectrum. We 
carry out numerical studies in the Glazek- Wilson and Wegner formalisms and 
present low lying eigenvalues and wavefunctions. We investigate the sensi- 
tivity of the spectra to various parameterizations of the similarity factor and 
other parameters of the effective Hamiltonian, especially the scale a. Our 
results illustrate the need for higher order calculations of the effective Hamil- 
tonian in the similarity renormalization scheme. 

PACS Numbers: ll.10.Ef, ll.10.Kk, 11. 10. St, 11.15.Tk, 12.38.Lg 



I. INTRODUCTION 



Any attempt to solve bound state problems in quantum field theory using Fock space 
based Hamiltonian methods immediately encounters a two-fold infinity problem. In a rela- 
tivistic system there are infinitely many energy scales and there are an infinite number of 
particles. Typical Hamiltonians of interest couple low energy scales with high energy scales 
which results in ultraviolet divergences. Furthermore, Hamiltonian couples every particle 
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number sector allowed by symmetries and at strong coupling brute force particle truncation 
can fail miserably. 

Light front Hamiltonian poses special problems for the treatment of ultraviolet diver- 
gences since the counterterms for transverse divergences in this case are in general expected 
to be nonlocal |I| in the longitudinal direction. In spite of the complexities due to renor- 
malization, light front approach to QCD is appealing from several physical aspects. First 
and foremost is the fact that in a cutoff theory where one imposes a cutoff on the mini- 
mum longitudinal momentum of the constituents, Fock vacuum is the exact eigenstate of 
the Hamiltonian and there is hope that one may find a direct link between QCD on one 
hand and constituent quark model and quark parton model on the other hand. 

In the best of all possible worlds, one can attempt to diagonalize the Hamiltonian in 
a single step. This is the spirit behind straightforward implementations of the Discretized 
Light Cone Quantization (DLCQ) method (for a review see Ref. ||). DLCQ has been quite 
successful in two dimensional models, but for QCD this approach may be quite ambitious. 
Recall efficient numerical algorithms designed for matrix diagonalization where the matrix is 
first brought to a tri-diagonal form which is then diagonalized [[|. In the same spirit, to keep 
the complexities under control, a two step process to the Hamiltonian diagonalization has 
been suggested The proposal is based on similarity renormalization scheme M devised 
independently by Glazek and Wilson and Wegner. 

In the first step, by integrating out the high energy degrees of freedom perturbatively 
one arrives at an effective Hamiltonian which is in a band diagonal form. At this step, 
one can identify the ultraviolet divergent part of the counterterms needed to be added 
to the Hamiltonian to remove ultraviolet divergences. It is quite advantageous to treat 
the ultraviolet divergences perturbatively especially in gauge theories since one can avoid 
pitfalls of other effective Hamiltonian approaches. For example, it is well-known that a 
simple truncation of the Fock space (like a Tamm-Dancoff truncation) leads to uncanceled 
divergences as a result of violations of gauge symmetry. 

In the second step, the effective Hamiltonian is diagonalized exactly. It is important to 
note that the high energy states integrated out contain small x gluons, i.e., gluons having 
small longitudinal momentum fraction. Since the vacuum is trivial, it is hoped that, as a 
result of integrating out small x gluons which are sensitive to long distance physics on the 
light front, the effective Hamiltonian may contain interactions responsible for low energy 
properties of QCD. Indeed Perry || found a logarithmic confining interaction in the qq 
sector in the lowest order effective interaction. 

Initial bound state studies in the similarity renormalization approach worked in either 
the non-relativistic limit |5| or in the heavy quark effective theory formalism |7|] to inves- 
tigate heavy-quark systems. Only recently, work has begun || in the context of glueball 
spectrum to address many practical problems, especially the numerical ones that one faces 
in this approach. Since the conceptual and technical problems one encounters in QCD are 
numerous, we have initiated a study of bound state problems in QCD in 2+1 dimensions 0. 
Our main motivation is not the fitting of data but a critical evaluation of the strengths and 
weaknesses of the various assumptions and approximations made in the similarity approach. 

In Ref. pi we studied the meson sector of 2+1 dimensional light front QCD using a 



Bloch effective Hamiltonian [10 in the first non-trivial order. The resulting two dimen- 



sional integral equation was converted into a matrix equation and solved numerically. The 
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vanishing energy denominator problem which leads to severe infrared divergences in 2+1 
dimensions was investigated in detail. We defined and studied numerically a reduced model 
which is relativistic, free from infrared divergences, and exhibits logarithmic confinement. 
The manifestation and violation of rotational symmetry as a function of the coupling was 
studied quantitatively. Our study indicated that in the context of Fock space based effective 
Hamiltonian methods to tackle gauge theories in 2+1 dimensions, approaches like similarity 
renormalization method may be mandatory due to uncanceled infrared divergences caused 
by the vanishing energy denominator problem. 

In this work we study similarity renormalization approach in the first non-trivial order to 
the same problem. The plan of this paper is as follows. In Sec. II we review the similarity 
renormalization approach in the Glazek- Wilson and Wegner formalisms. In Sec. Ill we 
present the effective Hamiltonian in the qq sector. Three different parameterizations of the 
similarity factor are discussed in Sec. IV. In Sec. V we present analytical calculations with 
the step function similarity factor which clearly shows the emergence of linear and square 
root confining interactions for large transverse and longitudinal separations respectively. In 
Sec. VI we perform numerical calculations and present low lying spectra and wavef unctions. 
Sec. VII contains discussion and conclusions. 



II. SIMILARITY RENORMALIZATION THEORY FOR THE EFFECTIVE 

HAMILTONIAN 

Realizing that a single step diagonalization of the Hamiltonian is too difficult, one may 
aim for the next best thing, namely, bringing the Hamiltonian to a band diagonal form. 
This is the method proposed by Glazek and Wilson and Wegner independently ||. 

Starting from a cutoff Hamiltonian Hb which includes canonical terms and counterterms 
we wish to arrive at an effective Hamiltonian H a defined at the scale a via a similarity 
transformation 

H a = S fj H B Si (2.1) 

where S a is chosen to be unitary. 

The boundary condition is Limit ._ K j H a = Hb- 

Introduce ant i- Hermit ian generator of infinitesimal changes of scale T a through 

S a = T ef™ da ' T °' (2.2) 

where T puts operators in order of increasing scale. 

For infinitesimal change of scale, S a = 1 — T a da and Si = 1 + T CT da. Then we arrive at 
the infinitesimal form of the transformation 

^=[H a ,T a }. (2.3) 

This equation which has been called the flow equation of the Hamiltonian is the starting 
point of the investigations. 

The basic goal of the transformation S a is that H a should be band diagonal relative to the 
scale a. Qualitatively this means that matrix elements of H a involving energy jumps much 
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larger than a should be zero. T a still remains arbitrary to a great extent. It is instructive 
to go through the steps of the derivation which leads to the Glazek- Wilson choice. 

We write H B = H BQ + H BI where H B0 is the free part and H BI is the interaction part of 
the bare cutoff Hamiltonian. A brute force way of achieving our goal is to define the matrix 
elements Hj ai j = f a ijH B iij where we have introduced the function f ai j = f(x a ij) with x a 
function of a 2 and AM|. The function f(x) should be chosen as follows: 

when a 2 » AM 2 , f(x) = 1 (near diagonal region); 

when a 2 << AM^, f(x) = (far off diagonal region); 

in between f(x) drops from 1 to (transition region). (2.4) 

Here AM|(= M 2 — M 2 ) denotes the difference of invariant masses of states i and j. Because 
of the properties of /, Hj ai j is band diagonal. What is wrong with such a choice of inserting 
form factors by hand at the interaction vertices? First of all, we have simply discarded 
degrees of freedom above a. Secondly, H a will have very strong dependence on a. Thirdly, 
to ensure that H a has no ultraviolet cutoff dependence, H B should contain canonical and 
counterterms. But, in light front Hamiltonian field theory, because of the complexities due 
to renormalization, a priori we do not know the structure of counterterms. 

Note that in the definition of H a given in Eq. ( |2.3|) the form of T a is still unspecified. In 
fact, a wide variety of choices are possible. In the following, we consider the choices made 
by Glazek and Wilson and Wegner. The price we have to pay for the use of flow equations 
is that it will generate complicated interactions even if the starting Hamiltonian has only 
simple interactions. For example, starting with a Hamiltonian which has only 2 particle 
interaction, the transformation will generate 3 particle interactions, 4 particle interactions, 
etc. 



A. Glazek- Wilson Formalism 



Writing H a = Hq + Hj a , noting that the free Hamiltonian Hq does not depend on a and 
taking matrix elements in free particle states, we have, 



[H a ,T a ]ij — (P i — Pj )T ai j + [Hi a ,T a ]ij 
where H \ i) — Pf \ i), etc. . i.e., 



1 dHj a ij 



[Hj a , TJ 



\P% Pj )P<rij- 



(2.5) 



(2.6) 



J aij OjO J crij J aij 

Since we want Hj a ij to be band diagonal, it is advantageous to trade for 
-f—Hj a ij which on integration has the chance to ensure that Hi a ij is band diagonal, 



_d_ 

da 

we use 



d 

da 



f 



-Hj 



v./ 



flii da 



1 dH Iai 
faij da 



(2.7) 



and arrive at 
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d 

da 



-H 



laij 



aij 



Pj )P<yij 



aij 



Jcnj J aij uu 



" /try • 



(2.8) 



Still T a ij is not denned. We next convert this equation into two equations, one denning the 
flow of Hi ai j and other defining T ai j. Recalling the starting equation Eq. ( |2.3| ) we add and 
subtract [H Ia) T a ]ij to the r.h.s. and arrive at 



d_ 

da 



-Hj 



v./ 



■ [H Ia ,T a ] i:j + — (Pr- P -)T aij 

Jcnj 

+ "7 (1 ~~ fcrij)[Hlcr,T a ]ij ~2 — 



(2.9) 



Glazek and Wilson choose T a to be 

1 



T, 



a.) 



pr - pr 

3 « 



d 



1 ferij) \Hlai P<r\ij , (hi faij)Hlaij 

da 



(2.10) 



Then from Eq. (|2.9| ), we have, 



dH 



laij 



da 



fcrij \Hla i Tq\ ij . 



Integrating Eq. (|2.11|) from a to oo, we arrive at, 



H 



laij 



f< 



/•OO 

HiBij — / da'[H Ia >,T a i 

J a 



'J 



(2.11) 



(2.12) 



Note that Hi^j is zero in the far off-diagonal region. This is clear from the solution 
given in Eq. (|2.12j ) since f(x) vanishes when x > 2/3. 

T ai j vanishes in the near diagonal region. When % is close to j, f ai j = 1 and both (l — f ai j) 
and 4z faij vanishes. It follows, then, from Eq. (|2.10|) that T ai j vanishes in the near-diagonal 
region. This guarantees that a perturbative solution to Hj ai j in terms of Hsnj will never 
involve vanishing energy denominators. 

Our next task is to derive the effective Hamiltonian to second order in perturbation 
theory. Using 



Hhlik — faikH B Iik 



and 



T, 



crkj 



PI ~ Pi 



d_ 

da 



(l n fcrkj)f<rkjHBIkj 



(2.13) 



(2.14) 



in Eq. ( [2.12|) , a straightforward calculation leads to 



Pfaij - H B RkHBIkj 



Qcrijk 



Qcrjik 



Pk-Pi 



Pk Pi 



(2.15) 
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where 



(2.16) 



We find that the effective Hamiltonian in similarity perturbation theory is a modification of 
the effective Hamiltonian in Bloch perturbation theory ||10| . 



B. Wegner Formalism 

In the Wegner formalism JTTJ] , the flow equation is given by 

^ = [r(l),H(l)]. (2-17) 

Wegner chooses 

r(l) = [H d ,H] = [H d ,H r ] (2.18) 

where Hj, is the diagonal part of the Hamiltonian and H r is the rest, i.e., H = H<i + H r . 
Here the word diagonal is used in the particle number conserving sense. It is important to 
note that is not the free part of the Hamiltonian and both and H r depend on the 
length scale 

The light front Hamiltonian has dimension of (mass) 2 and hence r has the dimension of 
(mass) 4 , I has dimension of -, — ^-tt. 

Expanding in powers of the coupling constant, 

H = H® + + Hf ] +H^ + ... (2.19) 

where the superscript denotes the order in the coupling constant, 

t(1) = [h£\hU] + [h£\hV} + ... . (2.20) 

Then, to second order, 

^ = [[H?\ H(%hP\ + [[HP,Hl% HW] + [[HP, Hl% tff ] +.... (2.21) 



Introduce |I2| the eigenstates of 

HP\i) = Pr\i). (2.22) 



Then, to second order 
dHuj 



ifi ~ PifH^ + H$% -(Pr- PrfH^ + . . . . (2.23) 



jt \- % ~ 3 > T%J ' I i ' r J*7 \ i 3 ' rl 3 

To first order in the coupling, 



6 



dK 



-\2rrW 



-{Pr-p-yH, 



which on integration yields 



(P7-P7V 

J. 3 



ri] ' 



(2.24) 



(2.25) 



where we have introduced the energy scale a via I = and used the fact that / = 
corresponds to the original bare cutoff. We notice the emergence of the similarity factor 



(p. -p. Y 



a i. j 



f, 

If we are interested only in particle number conserving (diagonal) part of the effective 
interaction, to second order we have, 



dH, 



dl 



Using 



(i) 



(Pr - Pr)H {1) 



the effective interaction generated to second order in the diagonal sector is 



H lij = H H ik H , 



{P- - P k ) + (pr - P k 



(pr-PD 2 +(p~-p k 



(2.26) 



(2.27) 



(2.28) 



Even though the second order formula is very similar to the one in Glazek- Wilson formalism 
when an exponential form is chosen for the similarity factor (see Sec. IV), we note a slight 
difference. In the Glazek- Wilson formalism, since the purpose is to bring the Hamiltonian 
into a band diagonal form, even in the particle number conserving sectors the large jumps in 
energies do not appear by construction. In the version of the Wegner formalism presented 
here the purpose is to bring the Hamiltonian in the block diagonal form in particle number 
sector so that large jumps in energies are allowed by the effective Hamiltonian. Note that 
small energy denominators do not appear in both formalisms. 



III. EFFECTIVE HAMILTONIAN IN THE QQ SECTOR 



In similarity renormalization approach due to Glazek and Wilson, to second order, the 
interacting part of the effective Hamiltonian in similarity renormalization approach is given 
by Eq. fl2.15|) . In this work we restrict ourselves to the qq sector. Then the states involved 
in the matrix elements i and j refer to qq states and k refer to qqg states. 

Following the steps similar to the ones outlined in Ref. 0, we arrive at the bound state 
equation 



M 



m 2 + k 2 



x(l — x) 



ip 2 (x, k)= SE ip 2 (x, k) 



C f dy 



2(2vr)< 



2(2vr) 2 f 
dq ^2{y,q) 



dy 
V 

KD' 



dq faij ip2{y,q) 



(x - y f 



(3.1) 
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Here x and y are the longitudinal momentum fractions and k and q are the relative transverse 
momenta. We introduce the cutoff rj such that r] < x,y < 1 — 77. We further introduce the 
regulator 5 such that | x — y \ > 5. Ultraviolet divergences are regulated by the introduction 
of the cutoff A on the relative transverse momenta k and q. The self energy contribution 



SE-- 



9' 



2(2tt) 2 ' 

-,2 



C f dy dq 6(x - y) [I- f aik ] xy 



2(k-q) 

(x-y) 



m 2 (x-y) 2 

x 2 y 2 



(ky — qx) 2 + m 2 (x — y) 2 
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2(2tt) 



dy dqd{y-x) [1 - f aik ] (1 - x)(l - y) 



x 



( q I k 1 2(q-k) \ 2 m 2 (y-x) 2 
{l-y "T 1-x "T" (y-x) ) "T (l-x)3(l-y)a 



[fc(l — y) — q(l — x)] 2 + m 2 (x — y) 1 
The boson exchange contribution 

V 9(x - y) 



ED 



(x-y) 



9crjik 
m 2 +q 2 _|_ (k-q) 2 

y {^-y) 



9uijk 



m 2 +k 2 
x 



x [K(k,x,q,y) + %Vi 



X) 



(y-x) 



9ojik 



+ 



m 2 +k 2 , (k-g) 2 
1—x x—y 



9crijk 



i-y 



2+fc 2 , (q-k) 2 q 2 +m 2 1 m 2 +q 2 , (q-k) 2 rr£ 
x (y-x) y l-y y-x 1 



+k 2 



X 



y,k,x) + iVj 



where 



K(k,x,q,y)~- 



'q | k 2 (k-q) \ . 
y x (x-y) J 
m 2 (x — y) 2 



+ 



k 2(k-q)' 



+ 



l-y 1-x (x-y) 



xy(l - x)(l - y)' 



(3.2) 



(3.3) 



(3.4) 



Tfl 



For all the / and g factors, 

2 _ 

x(l - x) " L ' y(l - y) 



, 2 k 2 + m 2 , ~ g 2 + m 2 

Mj = — and M = — -. (3.6) 



For x > y, M 2 k 



(k - qf 



m 



m 



x 



X 



(3.7) 



and 



for y > x, M k 



(q - kf 



+ 



q 2 + m 2 k 2 + m 2 



+ 



x 



X 



(3.8) 
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Before proceeding further, we perform the ultraviolet renormalization. The only ultra- 
violet divergent term arises from the factor 1 inside the square bracket in Eq. ( |3.2|) . We 
isolate the ultraviolet divergent term which is given by 



O ^divergent— ~ n , n N9 W / a V a Q / v 

2(2ir) 2 [Jo J -a xy(x — yy 

{2-x-yf 



f 1 r+ A (2 - 
+ / ay / dq - r— 

Jx+S J -A (1 — X)(l 

which is canceled by adding a counterterm. 



(1 - x)(l - y)(x - y)^ 



(3.9) 



IV. SIMILARITY FACTORS 



A. Parameterization I 



In recent numerical work JSJ , the following form for the similarity factor has been chosen: 



{AM 2 .) 2 



a i. j 



e 



(4.1) 



with AM?- = M? — M„ 2 where M 2 denotes the invariant mass of the state i, i.e., M 2 



. Then 



9auk= fan I da' f a , ik -Lu ik 



(AM 2 .) 2 



e 



(AM 2 ^ 2 



l jkJ 



(AMD* + (AM 2 ,) 2 



1 -e 



((AM 2 fc ) 2 + (AM 2 fc ) 2 



For the self energy contribution, i = j and we get 



1 

Qaijk Qajik Qaiik ~Z 



2(AM 2 fc ) 2 

1 -e ^~ 



(4.2) 



(4.3) 



Due to the sharp fall of / with a, the effective Hamiltonian has a strong dependence on a. 
Note that this parameterization emerges naturally in the Wegner formalism. 

B. Parameterization II 



St. Glazek has proposed the following form [13] for /, 



1 



(4.4) 



with 



9 



u. 



AM 2 



EM? + a 2! 



(4.5) 



uq a small parameter, and n„ an integer. The mass sum EM? = M? + M ? 2 . The derivative 



da 



2 n g 



2a 



EM? + a 2 \u 



u 



2 n 9 



2"s 



(l-tM,)*"* 1 


1 4 


( «<Tij(l-«0)\ 2 9 


2 






\«0(1— U ai j)) 





Note that for small u, both 1 — /(w) and j- vanish like w 2 ™ 9 . 



(4.6) 



C. Parameterization III 

For analytical calculations it is convenient to choose || a step function cutoff for the 
similarity factor: 



faij = 0(a 2 — AM?). 



(4.7) 



Then 



gaak = 0(° 2 ~ AM 2 ) fl(AM 2 - a 2 ) 6(AM 2 - AM 2 ). 



(4.8) 



It is the factor 8(AM 2 k — a 2 ) in g a ijk that prevents the energy denominator from becoming 
small. 



V. ANALYTICAL CALCULATIONS WITH THE STEP FUNCTION 

SIMILARITY FACTOR 

In this section we perform analytical calculations to understand the nature of the effective 
interactions generated by the similarity factor. Since there are no divergences associated with 
r) and A, we suppress their presence in the limits of integration in the following equations. 

A. Self energy contributions 



Consider the self energy contributions to the bound state equation Eq. 
the energy denominators to expose the most singular terms, we have, 



Rewriting 



SE-- 



9' 



2(2^ Wo dy 



2k? f L dy i dq 



0(x-6- y) 
x-y 

6(y-x-6) 
y-x 



{1 faik) 
{1 — faik} 



q . k _ 2(k-q) \ 1 , m 2 (x-y) 2 
y x (x—y) J x 2 y 2 



(fc-9) 2 1 q 2 +m 2 k 2 +m 2 

{x-y) y x 



\( <Z , k , 2(q-k)\ 


2 1 m 2 (y-x) 2 


[\l-y T l-x ~T (y-x) J 


-T (l_ a; )2 (1 _ y) 2j 



(q—k) 2 _|_ q 2 +rn 2 _ k 2 +m 2 



(5.1) 



y—x 



1-2/ 



l-x 
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The terms associated with 1 in the curly brackets lead to ultraviolet linear divergent terms 
which we cancel by counterterms (see Sec. III). They also lead to an infrared divergent term 
H which remains uncanceled. Explicitly this contribution is given by 

4 2(2^)2 ° f l dy J dq [ky -qx} 2 + m 2 (x-y) 2 

g 2 Ul 2 /"l r 1 

+4 2(2^ ° f L + 5 dy J dq [k(l-y)-q(l-x)} 2 + m 2 (x-y) 2 - (5 ' 2) 

This is simply indicative of the fact that terms associated with 1 in the curly bracket still 
has a vanishing energy denominator problem. We will address the resolution of this problem 
shortly. 

Let us next consider new infrared divergences that arise as a result of the modifications 
due to similarity factor. 



1. Leading singular terms 



Keeping only the most infrared singular terms in the numerators (i.e., for x > y, 4 ^_^2 

and for y > x, 4 ^~^ a ) and denominators (i.e., for x > y, -^E^tt an d for y > x ( q ~ k ^ 
have, 



SE 1 ~- 



fr 



2(2vr) 2 1 Jo 
c 2 

+ 



C f l dy I dq6(x-5-y) f aih 4 
9 



1 



(.V-b) 



we 



2 2tt 



(x - y) 2 

Cffcdy J dq6(y-x-5) f aik 4^ 



(y - x) 2 ' 



(5.3) 



The integral is given by 



dy I dq 



(x-y) 



:9 a 2 



{k - qf 
x-y 



+ / dy I dq 



x+5 



(y-x) 



:9 a 2 



(k-qf 
y-x 



(5.4) 



We change the transverse momentum variable, p = k — q. For x — 5 > y, we set x — y = z 
and for y > x + S we set y — x = z. Then, we have, 



: 2(2tt) 2 f 



x dz 



8 z z 



V 



1 ~ x dz 



^ dp9(**-^)+ -\d V B{c?-r- 



z^ 



p 



y/5 y/x Vl 



X 



(5.5) 



2. Sub-leading singular terms 
ingular terms 

similarity transformation. They are given by 



Next we study sub-leading singular terms containing in self energy generated by the 
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SE 2 = -4 



2(2tt) 2 



C f 



dy dq 9 a — 



1 



■k 2 q 2 



x-y 



(-- 



1 



Jx+8 J \ y — x J y — a; v 1 — x I — y 



x — y ^ x y ' (k — q) 2 
2 1 



yl (q-k) 2 



(5.6) 



where we have kept only (k — q) 2 term in the denominator since the rest vanish in the limit 
x — > y. As before, for x — 5 > y, we put x — y — z, k — q = p. With the symmetric 
integration in p, terms linear in p does not contribute. Only potential source of 5 divergence 
is the p 2 term in the integrand. Since p max = &\/z, after p integration - is converted into 
which is an integrable singularity. Same situation occurs for y > x. Thus there are no 
terms divergent in 5 coming from sub-leading singular terms. 



B. Gluon exchange contributions 

Let us next consider the effect of similarity factors on gluon exchange terms. 



1. Instantaneous gluon exchange 
From instantaneous interaction we have, 



V, 



inst 



-4 



9 



Cf j dy j dq tp 2 (y,q) fa 



2(2tt) 2 ~ j J ™» J ^ rzvy '" y Jm] {x- y y 
For the sake of clarity, it is convenient to rewrite this as 

Vi ns t = - 4 - 9 ^C f j dy J dq f aij tp 2 (y, q) 



!(2tt) : 
6{x-y-5) 



x— y v y x ' ^-y x> 



+ 



+ 



+ 



X -y { ik _ q) 2 + {3 ±_k± ){x _ y)+m 2 { l_l ){x _ y) 
x—y y t—y i—x' y t—y \—x> I 



(* - <?) 2 -(£-&)(*-!/)- m2 (ik - iX* - v) 

v y x ' y x y ' 



e(y-x-6) 



y-x 



y-x 



(q-k)' 



V 



- *)(y-x) + m 



2 ^-±)(y-*) 



y—x v 1— y \—x> v 1— y \—x> 



We have to seperately analyze the three types of terms in the numerator. 

First consider terms proportional to m 2 in the numerator. They are given by 

,2^2 



g tyi 1 I f 
- 4 ^y2 Cf J dy J dqf ^ My > q) 



X 



1 



+ 



1 



[ky — qx] 2 + m 2 (x — y) 2 [k(l — y) — q(l — x)] 2 + m 2 (x — y)' 



(5.7) 



(5.8) 



(5.9) 
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which leads to the logarithmic confining interaction in the nonrelativistic limit. Note, how- 
ever, that Eq. (|5.9|) is affected by a logarithmic infrared divergence arising from the vanishing 
energy denominator problem. The logarithmic infrared divergence is cancelled by the self 
energy contribution, Eq. ( |5.2|) . Thus we have explicitly shown that the logarithmically con- 
fining Coulomb interaction survives similarity transformation but the associated infrared 
divergence is canceled by self energy contribution. 

Next we look at the most singular term in the numerator in the limit x —>■ y. In this 
limit we keep only the leading term in the denominator and we get 



-4 



2(2vr)< 



C f dy dq ip 2 (y,q) f, 



0(x -y) 8{y - x) 
-y) 2 (y-x) 2 



(5.10) 



Lastly we look at the rest of the terms in the instantaneous exchange. Since we are interested 
only in the singularity structure, we keep only the leading term in the denominator and we 
get 



2 2kv 2Cf ! dy ! dqf ' 



2(2tt) 2 
(x-y) iq 2 (l-2y) 



x-y 



k 2 (l-2x) 



" n {k-q) 2 

9(y - x) vk 2 {\ - 2x) 



y(l-y) x(l-x) 



+ 



y-x 



q 2 {\ - 2y) 



x(l-x) y(l-y) 



(5.11) 



2. Transverse gluon exchange 
First, consider the most singular terms. 

Keeping only the most singular terms, the gluon exchange contribution is 



2(2tt) 
6(x-5- y) 



-Cf I dy I dq i[> 2 (y,q) x 



x-y 
9 2 



9ujik ~\~ 9crijk / A \i]^ 



(k-q) 2 



:-4) 



(x - yf 



+ 



2(2tt) 2 
(x-S-y) 



C f j dy f dq ^ 2 {y,q) x 

9{y - x - 5) 



9{y-x- S) 
y-x 



(q-W 
(y-x) 



(y 



X) 



(x - yf 



9crjik 9crijk 



+ 



(y - xf 



9crjik 9crijk 



Explicitly, for x > y, 



(5.12) 



9auk= 0{a 2 - M 2 ) 9{M% - M 2 k ) 9{M% - a 2 ), 
9^k= 9(a 2 - M 2 ) 9(M 2 k - M 2 ) 9(M 2 k - a 2 ). 



(5.13) 



We are interested in the situation x near y and i near j. Then 9(M 2 k — M 2 k ) 
9{M 2 k — M 2 k ) and 9(a 2 — M?-) = 1. Then the gluon exchange contribution is 
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A 2kr Cf l dy l dqMy,q) x 

- Qjx-S-y) f _J a 2_ (k-q) 2 
(x-y) 2 \ I 



x-y J J {y-x) { 



2 ( q -kr 

y-x 



(5.14) 



where we have used 6(x) = 1 — 8(—x). Combining with the most singular part of the 
instantaneous contribution given in Eq. ( |5.10| ) we arrive at 



2(2vr)< 



C f dy dq ^ 2 {y,q) x 



6(x-5-y) 
(x - y) 2 



9 a 2 - 



(k-q) 2 \ 9(y-x-S) 



x-y 



(y - xf 



9 a 2 



(g ~ kf 
y-x 



(5.15) 



For convenience we change variables. For x > y, we put x — y = p+ and k — q = p 1 and for 
y > x, we put y — x = ^+ and q — k = p 1 where P + is the total longitudinal momentum. 
Thus we arrive at 



-4 



2(2tt) 



C f P^ 



r r p+x p + if 

Jdpj p+s dp + ^ 2 {x-—,k-p 1 )-^9ya 2 



{p 1 ) 2 P + * 



P 



+ / dp 



P+8 

Consider the Fourier transform 

,2 



P+(l-x) n + 1 / (n 1 ) 2 P + v 



P^ 



(P 



P 



(5.16) 



V(x-,x^)= -4 







2(2vr 



P+5 (p 



P+5 



P max 



P 



dp 1 e* p x ~ ip x 9 



a' 



{p 1 ) 2 P + ' 



P 



(5.17) 



where p x max = cry We are interested in the behavior of V(x , x 1 ) for large x ,x l . For 
large x~ , nonnegligible contribution to the integral comes from the region q + < j^py. For 

large x 1 , we need p^x 1 to be small, i.e., {p l max ) 2 < t^tw, i.e., p + < {x i^ a 2 ■ Thus we have 



the requirements, p + < npr, p + < ^ 2(7 2 ■ We make the approximations 



dp e 



1 — ip x 



2P 1 

Prn 



(5.18) 



and e2 q+x « 1. 

For large we have p + < rpr < , ^3 g2 , the upper limit of p + integral is cut off by ^py. 
Adding the contributions from both the integrals (which are equal), for large x~ , we have 

.2 j 



2(2tt) 



P+ x- = 



(5.19) 



Thus for large a; the similarity factors have produced a square root potential but it is also 
infrared singular. 
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For large x l the upper limit of p + integral is cut off by -n^w^i and we get, 



Vi^x-.x 1 ) 



32 



C f a 



x a 



1 

71 



(5.20) 



2(2tt) 2 

For large x 1 , similarity factors have produced a linear confining potential which is also 
infrared singular. We note that the rotational symmetry is violated in the finite part of the 
potential. In both cases, however, the infrared singular part is ~32 2 ^ 2 Cf o-^ which is 
exactly canceled by the infrared contribution generated by similarity transformation from 
self energy, Eq. (157 



Lastly, we consider the terms that go like Keeping only the leading term in the 

x y 

energy denominator, we have, 



r 



C f dy dq 4> 2 (y,q) x 



2(2tt) 2 

0(x-y) g<rjik + 9aijk[k 2 (l-2x) q 2 {l-2y) 



x-y 



(k - qf 



x(l — x) 



2/(1-2/) 



+ 



6(y-x) g ajik + g a ijk\q 2 (l-2y) k 2 (l-2x) 



2/(1 -y) 



x(l — x) 



y — x (q — k) 2 
With the step function cut off we have 

Qajik + g*ijk ~ f m jO(AM 2 k - a 2 ). 
Then, combining Eq. (|5.11|) and Eq. Q5.21]) for the sub-leading divergences, we get, 



(5.21) 



(5.22) 



-2 



2(2tt) 



C f dy dq f aij tp 2 (y,q) x 



0(x - y) 1 
x — y (k — q) 2 
6{y - x) 1 







a 



(k-q) 2 \rk 2 (l-2x) q 2 (l-2y) 



x-y 



y — x (q — k)'- 



x(l-x) y(l-y) 

(q-k) 2 \ rq 2 (l-2y) k 2 (l-2x) 



2/(1 -y) x(l-x) 



(5.23) 



y-x 

Taking the Fourier transform of this interaction, a straightforward calculation shows that 
no log 5 divergence arise from this term. 



3. Summary of divergence analysis 

The logarithmic confining Coulomb interaction of 2+1 dimensions is unaffected by sim- 
ilarity transformation and is still affected by a logarithmic divergence which is however 
cancelled by a logarithmic divergence from self energy contribution. Similarity transforma- 
tion leads to a non-cancellation of the most singular ( / x _ y \i ) term between instantaneous 
and transverse gluon interaction terms. This leads to a linear confining interaction for large 
transverse seperations and a square root confining interaction for large longitudinal sepa- 
rations. However, non-cancellation also leads to 4j divergences where 5 is the cutoff on 
| x — y | . This divergence is cancelled by new contributions from self energy generated by 
similarity transformation. The subleading singular terms do not lead to any divergence 

x y 

in 5. 
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VI. NUMERICAL STUDIES 



The integral equation is converted in to a matrix equation using Gaussian Quadrature. 
The matrix is numerically diagonalized using standard LAPACK routines |14j. (For details 
of numerical procedure see Ref. f9j.) With the exponential form and the step function form 
of the similarity factor, the integral over the scale in the definition of g factors in Eq. (|2.16|) 
can be performed analytically as shown in Sec. IV. For parametrization II, we perform the 
integration numerically using ns quadrature points. 

The first question we address is the cancellation of divergences which are of two types: 
(1) the In 5 divergence in the self energy and Coulomb interaction which has its source 
in the vanishing energy denominator problem that survives the similarity transformation 
and (2) A= divergences in the self energy and gluon exchange generated by the similarity 
transformation. In Table I we present the 5 independence of the first five eigenvalues of the 
Hamiltonian for g=0.6. Results are presented for three parametrizations of the similarity 
factor, namely, the exponential form, the form proposed by St. Glazek and the step function 
form used in our analytical studies. It is clear that the Gaussian Quadrature effectively 
achieves the cancellation of 5 divergences. Recall that in the study || of the same problem 
using Bloch approach, negative eigenvalues appeared for g=.6 when 5 was sufficiently small 
(for example, .001) which was caused by the vanishing energy denominator problem. Our 
results in the similarity approach for the same coupling shows that this problem is absent 
in the latter approach. 

Next we study the convergence of eigenvalues with quadrature points. In Table II we 
present the results for all three parametrizations of the similarity factor for g=0.2 with the 
transverse space discretized using k = ^tan^- where w's are the quadrature points |J. The 
table show that for m — 1, convergence is rather slow for all three choices of the similarity 
factor compared to the results in Bloch approach. Among the three choices, parameterization 
II shows better convergence. 

Let us next discuss the nature of low lying levels and wavefunctions. First we show the 
ground state wavefunctions for all three similarity factors for a given choice of parameters 
in Fig. 1. As is anticipated, step function choice produces a non-smooth wavefunction. For 
parameterizations I and II, the wavefunctions show some structure near x = 0.5. From our 
previous experience with calculations in the Bloch formalism, we believe that the structures 
indicate poor convergence with the number of grid points. 

Consider now the structure of low lying levels. Recall that in the Bloch formalism, 
the ordering of levels was I — 0, 1, 0, ... corresponding to logarithmic potential in the 
nonrelativistic limit. In the presence of effective interactions generated by the similarity 
transformation, obviously the level ordering changes. Now we have additional confining 
interactions which, however, act differently in longitudinal and transverse directions. The 
wavefunctions are presented in Fig. 2a for parametrization II for the first four low lying 
levels. 

There is extra freedom in parametrization II due to the presence of in the def- 

inition of u a ij, Eq. ( [4.5| ). For zero transverse momentum of constituents, SM^ has the 
minimum value 8m 2 . Thus relative insensitivity of parametrization to a in parametrization 
II for small values of a may be due to this factor. When we consider the heavy fermion 
mass limit, presence of 8m 2 in u a ij enhances the effect of similarity factor. In Fig. 2b we 
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present the wavefunctions corresponding to first four levels for parametrization II with 8m 2 
subtracted from EM^ in u a ij. From Figs. 2a and 2b note that the fourth level is different 
for parametrization II with and without 8m 2 in u a ij. 

By suitable choice of parameters we can study the interplay of rotationally symmet- 
ric logarithmically confining interaction and effective interactions generated by similarity 
transformation. Since for a given g, strength of the logarithmic interaction and similarity 
generated interactions are determined by m and a respectively, for m » a we should re- 
cover the Bloch spectrum ||. Upto what levels the recovery occurs, of course depends on 
the exact value of m. As we already observed, for parametrization II this will happen only if 
8m 2 is subtracted from EM 2 •. For this case, we present the first four levels for m — 10 and 
a = 4 in Fig. 3 which clearly shows the level spacing corresponding to the Bloch spectrum. 

Finally, we discuss the sensitivity of the spectra to the similarity scale a. Ideally the 
low lying energy levels should be insensitive to a. However we have calculated the effective 
Hamiltonian to only order g 2 and we expect significant sensitivity to a. In Tables III and IV 
we present the lowest five eigenvalues for all three parametrizations of the similarity factor 
for g = 0.2 and g = 0.6 respectively. As expected a dependence is greater for larger value 
of g. Among the three parametrizations, the paramterization II is least sensitive to a. In 
order to check whether this behaviour is due to the presence of 8m 2 in EM 2 - in the definition 
of u a ij we present the results in Table V for parametrization II with 8m 2 subtracted from 
EM 2 -. It is clear that sensitivity to a is still considerably less compared to the other two 
parametrizations. This may be due to the fact that EM 2 - is added to a 2 in the definition of a. 
Note that in parametrization II sensitivity to a is controlled also by additional parameters 
uo and n g . The sensitivity to parameters n g and uq are presented in Tables VI and VII 
respectively. 

VII. SUMMARY, DISCUSSION AND CONCLUSIONS 

An attempt to solve 2+1 dimensional gauge theories using the Bloch effective Hamilto- 
nian has revealed |J problems due to uncancelled infrared divergences. They arise out of 
vanishing energy denominators. Similarity renormalization formalism attempts to solve the 
bound state problem in a two step process. At the first step, high energy degrees of freedom 
are integrated out and ultraviolet renormalization carried out perturbatively. In the second 
step, the effective Hamiltonian is diagonalized non-perturbatively. In this work, we have 
studied the bound state problem in 2+1 dimensional gauge theories using the similarity 
approach. 

In order to have a better understanding of the numerical results, we have performed 
analytical calculations with step function form for the similarity factor. Many interesting 
results emerge from our analytical calculations. First of all, it is shown that due to the 
presence of instantaneous interactions in gauge theories on the light front, the logarithmic 
infrared divergence that appeared in the Bloch formalism persists in two places, namely 
a part of the self energy contribution and the Coloumb interaction that gives rise to the 
logarithmically confining potential in the nonrelativistic limit. However the terms that 
persist are precisely those that produce a cancellation of resulting infrared divergences in the 
bound state equation. The rest of the infrared problem that appeared in the Bloch formalism 
due to the vanishing energy denomnator problem is absent in the similarity formalism. 
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Similarity transformation however prevents the cancellation of the most severe ^_^ 2 
singularity between instantaneous gluon exchange and transverse gluon exchange interac- 
tions and produces divergences in the self energy and gluon exchange contributions which 
cancells between the two in the bound state equation. The resulting effective interaction 
between the quark and antiquark grows linearly with large transverse separation but grows 
only with the square root of the longitudinal separation. This produces violations of rota- 
tional symmetry in the bound state spectrum. We have also verified that no In 5 divergence 
result from the —7- singularity in the self energy and gluon exchange contributions. 

x y 

In the Glazek- Wilson formalism the exact form of the similarity factor f a is left unspec- 
ified. In the literature an exponential form has been used in numerical calculations ||. For 
analytical calculations it is convenient to choose a step function even though it is well-known 
that it is not suitable for quantitative calculations [15| . There is also a proposal due to Stan 
Glazek which has two extra free parameters. We have tested all three parametrizations in 
our work. Our numerical results indeed show that step function choice always produces non- 
smooth wavef unctions. Parameterization II costs us an extra integration to be performed 
numerically but convergence is slightly better for small g compared to exponential form. All 
three parametrizations produce violations of rotational symmetry even for small g. When an 
exponential form is used in the Glazek- Wilson formalism, the resulting effective Hamiltonian 
differs from the Wegner form only by an overall factor that restricts large energy diffrences 
between initial and final states. Numerically we have found this factor to be insignificant. 

We have studied the sensitivity of the low lying eigenvalues to the similarity scale a. 
Since the effective Hamiltonian is calculated only to order g 2 results do show sensitivity to 
o. Among the three parametrizations the form II is least sensitive to a due to the functional 
form chosen. We have also studied the sensitivity of eigenvalues to the parameters u and 
n g . 

The bound state equation has three parameters m, g 2 and a with dimension of mass. The 
strength of the logarithmically confining interaction is determined by m and the strength of 
the rotational symmetry violating effective interactions generated by similarity transforma- 
tion is determined by a. For a given g we expect the former to dominate over the latter for 
m >> a. An examination of low lying eigenvalues and corresponding wavefunctions show 
that this is borne out by our numerical calculations. 

A major problem in the calculations is the slow convergence. Compared to the Bloch 
formalism, in calculations with the similarity formalism various factors may contribute to this 
problem with the Gauss quadrature points. One important factor is the presence of linear 
confining interactions generated by the similarity transformation. It is well known that such 
interactions are highly singular in momentum space. Another factor is the presence of ^= 
divergences the cancellation of which is achieved numerically. It is of interest to carry out 
the same calculations with numerical procedures other than the Gauss quadrature. However, 
one should note that calculations in 3+1 dimensions employing basis functions and splines 
have also yielded |§ wavefunctions which show non-smooth structures. 

An undesirable result of the similarity transformation carried out in perturbation theory 
is the violation of rotational symmetry. Our results show that this violation persists at all 
values of g for m = 1. Such a violation was also observed in 3+1 dimensions. In that 
case the functional form of the logarithmic potential generated by similarity transformation 
is the same in longitudinal and transverse directions but the coefficients differ by a factor 
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of two. Same mechanism in 2+1 dimensions makes even the functional forms different. 
The important questions are whether the confining interactions generated by the similarity 
transformation are an artifact of the lowest order approximation and if they are not, then, 
whether the violation of rotational symmetry will diminish with higher order corrections to 
the effective Hamiltonian. Recall that the high energy degrees of freedom has been integrated 
out and the effective low energy Hamiltonian determined only to order g 2 . A clear answer 
will emerge only after the determination of the effective Hamiltonian to fourth order in the 
coupling. 
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FIGURES 



FIG. 1. The ground state wavefunction for different choices of the similarity factor using 
the tan parametrization for transverse momentum grid and for n\ = 40, n 2 = 80, m—1.0, 
g = 0.2, a = 4.0 and r] — 8 — 10~ 5 , re = 10.0 as a function of x and k. (A) Parametrization 
I, (B) Parametrization II with n g — 2,u — 0.1, ns = 500, (C) Parametrization III. 



FIG. 2a. The wavefunctions corresponding to the lowest four eigenvalues as a function of x 
and k for parametrization II. The parameters are as in FIG. 1. (A) Lowest state, (B) first 
excited state, (C) second excited state, (D) third excited state. 



FIG. 2b. Same as in FIG. 2a but with 8m 2 subracted from EM, 2 ,- . 



FIG. 3. The wavefunctions corresponding to the lowest four eigenvalues as a function of 
x and k with parametrization III with 8m 2 subtracted from EM 2 -. The parameters are 
g = 0.2, 7] = 5 = .00001, m = 10.0, re = 10.0, n g = 2, u = 0.1, m = 40, n 2 = 80. (A) 
Lowest state, (B) first excited state, (C) second excited state, (D) third excited state. 
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r 



Parametrization I 


0.1 


4.89535 


4.90359 


4.90420 


4.90420 


4.90482 


0.01 


5.62612 


6.38083 


6.82963 


7.20037 


7.36414 


0.001 


5.68417 


6.42147 


6.90879 


7.30609 


7.64650 


0.0001 


5.68432 


6.42148 


6.90909 


7.30611 


7.64677 


0.00001 


5.68432 


6.42148 


6.90909 


7.30611 


7.64677 





Parametrization II 


0.1 


4.55364 


4.55668 


4.55668 


4.55668 


4.55669 


0.01 


4.86066 


5.33491 


5.49693 


5.59838 


5.79111 


0.001 


4.87607 


5.35671 


5.59226 


5.64476 


5.88613 


0.0001 


4.87604 


5.35671 


5.59236 


5.64477 


5.88615 


0.00001 


4.87604 


5.35671 


5.59236 


5.64477 


5.88615 





Parametrization III 


n i 

U. 1 


5.12410 


5.13039 


5.13101 


5.13101 


5.13754 


0.01 


6.02600 


6.94326 


7.50445 


7.98054 


8.39927 


0.001 


6.00968 


6.97160 


7.55376 


8.07749 


8.49199 


0.0001 


5.96636 


6.97160 


7.51814 


8.07751 


8.46524 


0.00001 


5.96636 


6.97160 


7.51814 


8.07751 


8.46524 



TABLE I. Variation with 5 of the first five eigenvalues of the full hamiltonian (excluding 
the less significant imaginary term). The parameters are m = 1.0, g = 0.6, n\ = 58, 
n 2 = 58, 7] = 0.00001, A = 20.0, a = 4.0, (u = 0.1, n g = 2 and ns = 500 (for a integration) 
in parametrization II) 
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Tlx 


n-2 


Parametrization I 


10 


10 


4.320 


4.353 


A O C7 

4.357 


A O C7 

4.357 


A OP 1 

4.361 


20 


20 


4.375 


4.442 


4.484 


4.484 


4.485 


20 


30 


4.398 


4.482 


4.546 


4.583 


4.610 


20 


40 


4.411 


4.502 


4.570 


4.615 


4.656 


30 


40 


4.412 


4.503 


4.570 


4.615 


4.655 


40 


50 


4.420 


4.515 


4.585 


4.634 


4.678 


40 


60 


4.426 


4.524 


4.594 


4.645 


4.692 


40 


80 


4.434 


4.535 


4.607 


4.661 


4.709 


Tli 




Parametrization II 


10 


10 


4.163 


4.194 


4.194 


4.194 


4.203 


20 


20 


4.186 


4.244 


4.276 


4.276 


4.277 


20 


30 


4.192 


4.256 


4.296 


4.323 


4.329 


20 


40 


4.195 


4.262 


4.304 


4.335 


4.344 


30 


40 


4.195 


4.262 


4.304 


4.335 


4.344 


40 


50 


4.197 


4.266 


4.308 


4.341 


4.353 


40 


60 


4.199 


4.268 


4.311 


4.345 


4.359 


40 


80 


4.201 


4.272 


4.315 


4.350 


4.367 


Hi 




Parametrization III 


10 


10 


4.360 


4.379 


4.381 


4.381 


4.391 


20 


20 


4.469 


4.533 


4.572 


4.572 


4.572 


20 


30 


4.503 


4.604 


4.674 


4.721 


4.749 


on 


a n 


4.520 


4.632 


4.703 


4.768 


4.810 


30 


40 


4.528 


4.636 


4.714 


4.768 


4.811 


40 


50 


4.542 


4.657 


4.736 


4.797 


4.848 


40 


60 


4.548 


4.668 


4.748 


4.813 


4.866 


40 


80 


4.556 


4.683 


4.764 


4.833 


4.889 



TABLE II. Convergence of eigenvalues with rii and ri2 for the parametrization k = 
Han(qn/2). The parameters are m=1.0, g=0.2, r] = 0.00001, 5 = 0.00001, re = 10.0, 
a = 4.0, (-u =0.1, n g =2 and the number of quadrature points ns = 500 for a integration for 
parametrization II). 
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Eigenvalues (M 2 ) 


a 


Parametrization I 


2.0 


4.214 


4.285 


4.329 


4.365 


4.393 


4.0 


4.434 


4.535 


4.607 


4.661 


4.709 


6.0 


4.701 


4.821 


4.914 


4.980 


5.043 


a 


Parametrization II 


2.0 


4.157 


4.214 


4.248 


4.260 


4.276 


4.0 


4.201 


4.272 


4.315 


4.350 


4.367 


6.0 


4.266 


4.350 


4.404 


4.446 


4.483 


a 


Parametrization III 


2.0 


4.254 


4.346 


4.395 


4.443 


4.477 


4.0 


4.556 


4.683 


4.764 


4.833 


4.889 


6.0 


4.927 


5.075 


5.179 


5.262 


5.333 



TABLE III. Variation with a of the full hamiltonian (excluding the imaginary term). 
The parameters are m = 1.0, g = 0.2, nl = 40, n2 = 80, rj = 0.00001, k = 10.0, 
(k = Han(q7r/2)), 5 = 0.00001, (u = 0.1, n g = 2, and ns = 500 for a integration for 
parameterization II) 
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Eigenvalues (M 2 ) 


a 


Parametrization I 


2.0 


4.941 


5.418 


5.708 


5.713 


5.950 


4.0 


5.684 


6.421 


6.909 


7.306 


7.647 


6.0 


6.678 


7.617 


8.274 


8.801 


9.267 


a 


Parametrization II 


2.0 


4.769 


5.152 


5.226 


5.372 


5.485 


4.0 


4.876 


5.357 


5.592 


5.645 


5.886 


6.0 


5.071 


5.655 


6.016 


6.159 


6.316 


a 


Parametrization III 


2.0 


4.888 


5.589 


5.882 


6.120 


6.263 


4.0 


5.966 


6.972 


7.518 


8.077 


8.465 


6.0 


7.359 


8.603 


9.360 


10.083 


10.621 



TABLE IV. Variation with a of the full hamiltonian (excluding the imaginary term). The 
parameters are m = 1.0, g = 0.6, nl = 58, n2 = 58, r] = 0.00001, A = 20.0, (k = (1 _^ +m ), 
5 = 0.00001, (uq = 0.1, n g = 2, and ns = 500 for a integration for parameterization II) 
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Eigenvalues (M 2 ) 


9 


a 


Parametrization II 




2 


4.126 


4.166 


4.176 


4.189 


4.201 


0.2 


4.0 


4.172 


4.235 


4.273 


4.298 


4.304 




6.0 


4.241 


4.321 


4.371 


4.411 


4.445 




2.0 


4.739 


5.003 


5.020 


5.134 


5.193 


0.6 


4.0 


4.802 


5.222 


5.350 


5.470 


5.626 




6.0 


4.993 


5.541 


5.876 


5.940 


6.156 



TABLE V. Variation with a of the full hamiltonian (excluding the imaginary term) 
after subtracting 8m 2 from SM 2 - in the definition of u ai j. The parameters are m = 1.0, 
rj = 0.00001, 5 = 0.00001, u = 0.1, n g = 2, and ns = 500 

1) for g = 0.2, k = Han(qn/2) with k = 10.0 and nl = 40, n2 = 80 

2) for g = 0.6, k = ^_ff +m with A = 20.0, and nl = n2 = 58. 
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9 


a 


M 2 (Parameterization II) 




2.0 


4.781 


5.049 


5.051 


5.161 


5.221 


0.6 


4.0 


4.843 


5.235 


5.391 


5.464 


5.639 




6.0 


5.030 


5.538 


5.848 


5.956 


6.106 



TABLE VI. Variation with a of the full hamiltonian (excluding the imaginary term) after 
subtracting the 8m 2 term from EM 2 - in the definition of u a ij. The parameters are m = 1.0, 

v = 0.00001, 5 = 0.00001, u = 0.1, n g = 1, and ns = 500 g = 0.6, k = (1 _ffig+m with 
A = 20.0, and nl = n2 = 58. 
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Uq 


M 2 (Parameterization II) 




0.2 


4.990 


5.548 


5.887 


5.969 


6.171 




0.1 


4.802 


5.222 


5.350 


5.469 


5.626 


0.6 


0.05 


4.746 


5.076 


5.091 


5.253 


5.313 




0.01 


4.813 


5.003 


5.029 


5.069 


5.130 



TABLE VII. Variation with uq of the full hamiltonian (excluding the imaginary term) after 
subtracting the 8m 2 term from SM 2 - the definition of u ai j. The parameters are m = 1.0, 
rj = 0.00001, 5 = 0.00001, a = 4.0, n g =2, and ns = 500 g = 0.6, k = (1 _^" +m with 
A = 20.0, and nl = n2 = 58. 
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Figure 2b 




30 




31 



